*** Import 

use "$data/figure_A1_input.dta", clear

		twoway (scatter vdsa_hyv_share hyv_share, mcolor(black) msymbol(circle_hollow) ytitle(Share of area under HYV ( VDSA ) ) ///
		xtitle(Share of area under HYV ( AIS )) msize(medsmall)) || (lfit vdsa_hyv_share hyv_share, lcolor(blue) legend(off)) , ylabel(,nogrid) graphregion(color(white)) || (function y = x, ra(hyv_share) clpat(dash) lcolor(gray)) 
		graph export "$results/graphs/figure_A1.pdf", as(pdf) replace
		reg vdsa_hyv_share hyv_share, robust
		sum vdsa_hyv_share hyv_share if hyv_share !=., d

